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Orbital and Mass Ratio Evolution of Protobinaries Driven by 

Magnetic Braking 

Bo Zhao 1 , Zhi-Yun Li 1 
ABSTRACT 

The majority of stars reside in multiple systems, especially binaries. The 
formation and early evolution of binaries is a longstanding problem in star for- 
mation that is not yet fully understood. In particular, how the magnetic field 
observed in star-forming cores shapes the binary characteristics remains relatively 
unexplored. We demonstrate numerically, using an MHD version of the ENZO 
AMR hydro code, that a magnetic field of the observed strength can drastically 
change two of the basic quantities that characterize a binary system: the orbital 
separation and mass ratio of the two components. Our calculations focus on 
the protostellar mass accretion phase, after a pair of stellar "seeds" have already 
formed. We find that, in dense cores magnetized to a realistic level, the angular 
momentum of the material accreted by the protobinary is greatly reduced by 
magnetic braking. Accretion of strongly braked material shrinks the protobinary 
separation by a large factor compared to the non-magnetic case. The magnetic 
braking also changes the evolution of the mass ratio of unequal-mass protobinaries 
by producing material of low specific angular momentum that accretes preferen- 
tially onto the more massive primary star rather than the secondary. This is 
in contrast with the preferential mass accretion onto the secondary previously 
found numerically for protobinaries accreting from an unmagnetized envelope, 
which tends to drive the mass ratio towards unity. In addition, the magnetic 
field greatly modifies the morphology and dynamics of the protobinary accre- 
tion flow. It suppresses the traditional circumstellar and circumbinary disks 
that feed the protobinary in the non-magnetic case; the binary is fed instead by 
a fast collapsing pseudodisk whose rotation is strongly braked. The magnetic 
braking-driven inward migration of binaries from their birth locations may be 
constrained by high-resolution observations of the orbital distribution of deeply 
embedded protobinaries, especially with ALMA and JVLA. 

Subject headings: binary, accretion disks — magnetic fields — stars: formation 
— magnetohydrodynamics (MHD) 
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Introduction 



The majority of the Galactic field stars reside in multiple systems (|Duauennov fc Mayor 



1991 



2005 



Fischer fc Marcv 



Raghavan et al. 



(jDuquennoy Sz Mayor 



1992 



2010 



Mason et al 



Janson et al 



1991 



1998 



Shatsky fc Tokovininll2002t iGoodwin fc Kroupa 
201211; most of those systems (~ 75%) are binaries 



Tokovinin fc Smekhqv 2002h. The multip l icity of young stellar 



objects is even highe r ( Reipurth fc Zinnecker 19931 ; Mathieu et al. 2000 ; Duchene et al. 2004 
Duchene et al.l 120071 ). indicating that the formation of multiple systems, especially binaries, 
is a major, perhaps the dominant, mode of star formation. 

How binaries form and evolve remains uncertain. This is particularly true for the ear- 
liest, Class 0, phase, when the stellar seeds are still deeply embe dded inside a massi ve en- 
velope. Although millimeter interferometric observations such as lLooney et al.l ( 120001 ) have 
uncovered some Class binaries, their number is still small, and the distr ibution of orbital 
separation of such protobinaries is not well constrained. iMaury et al.l ( 120101 ) surveyed 5 Class 
sources with PdBI at sub-arcsecond resolution, and found a surprising result: only one of 
the sources has a potential protostellar companion, and t he companion is ~ 1 , 900 AU away 



from t he primary. Combining their sample with that of lLooney et al.l ( 120001 ). IMaury et al. 



(120101 ) found that none of the 14 Class sources in the combined sample has a companion 
with separation between ~ 150 — 550 AU, which is incon sistent with the binary fraction of 
~ 18% for Class I sources in the same separation range (jConnelley et al.ll2008l ). A similar 
difference was found i n the CARMA sur vey of 6 Class and 3 Class I sources in the Ser- 
pens molecular cloud (jEnoch et al.ll201ll ). None of the Class objects has any detectable 
protostellar companion closer than 2,000 AU (down to the resolution limit ~ 250 AU), 
whereas one Class I object has a companion at a projected distance of ~ 870 AU. Although 
the statistical significance of the difference is still relatively low, it brings into sharp focus 
the possibility that the distribution of orbital separation of protobinaries may differ signifi- 
cantly from that of their more mature counterparts. In other words, binaries may migrate 
substantially from their birth locations, a situation somewhat analogous to that inferred for 
hot Jupiters, although it is unclear whether protobinaries would migrate inward or outward 
to fill the apparent gap between ~ 150 — 550 AU. The direction (inward vs outward) of pro- 
tobinary orbital migration should be better constrained observationally with large surveys 
of deeply embedded sources using ALMA and JVLA. 

This paper focuses on one possible mechanism for moving the protobinaries away from 
their birth locations: magnetic braking. It may seem counter-intuitive that magnetic fields 
can change the binary orbit, because the magnetic forces do not act on the stars directly. 
However, the orbits of protobinaries are determined mainly by the angular momentum of 
the material to be accreted, which can be strongly affected, perhaps even controlled, by 
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the magnetic field, through magnetic braking. The main goal of this paper is to quantify 
the extent to which the protobinary orbit is modified by a magnetic field of the observed 
strength. 

The strength of magnetic fields in star formation is usually measured by the dimension- 
less mass-to-flux ratio A. It is the mass of a regi on divided by the magnetic flux threading the 
region in units of the critical value (2ir G 1 / 2 )~ 1 (Nakamura Sz Nakanolll978l ). For a sample of 
dense cores in the nearby dark clouds, iTroland fc Crutcherl (120081 ) inferred a mean value for 
\i os ~ 4.8, based on the field strength and column density along the line-of-sight. Geomet- 
ric corrections for projection effects should reduce this value by a factor of 2-3, yielding an 
intrinsic value of A of a few typically. Such a magnetic field is generally not strong enough 
to prevent the dense core from gravitational collapse and star formation. It is, however, 
strong enough to affect the angular momentum evolution of the collapsing core in general 
and binary formation in particular. 

There have been several studies of the magnetic effects on binary formation (e.g., 
Price fc Batdl2007l ; iHennebelle fc Teyssierll2008l ; iMachida et al.ll2010l ). These studies focused 
primarily on the classical mode of binary formation inside an isolated rotating core, through 



core frag mentation ind uced by an m = 2 density perturbation (e.g., iBoss fc Bodenheimer 
19791 . see iKratterl |20 1 ll for a recent review). A general conclusion is that the fragmentation 



can be supp ressed by a rather weak mag netic field if the core is only weakly perturbed. 
For example, IHennebelle fc Teyssierll2008l found that, for a density perturbation of 10%, the 
fragmentation is suppressed by a magnetic field corresponding to A = 20 (see their Fig. 3), 
much weaker than the observationally inferred field. In dense cores that are magnetized to a 
more realistic level (with A ~ a few), fragmenta t ion can still occur, although a large ampli- 
tude perturbation is needed (IPrice &: Bate! 120071 ; IHennebelle &: Teyssierl 120081 ) . A limitation 
of the existing MHD calculations is that they are confined mostly to the pre-stellar phase of 
core collapse and fragmentation leading up to the formation of two binary seeds. How the 
binary seeds evolve during the subsequent protostellar mass accretion phase in the presence 
of a dynamically important magnetic field remains little explored. It is the focus of our 
investigation. 

Our investigation of the magnetic effects on protobinary evolution will be carried out 
using an MHD version o f the ENZO AMR hydro code that includes a sink particle t reatment 
farvan fc Normanlll997l : b'Shea et allbooJ IWang fc Abell2009l : IWang et alJboioh . In § El 
we discuss the setup for the initial binary seeds and the rotating magnetized protobinary 
envelope to be accreted by the seeds, as well as the numerical code used. In § [31 we present 
numerical results for the evolution of initially equal-mass binary seeds. We find that a 
magnetic field of the observed strength can shrink the protostellar orbit by more than an 
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order of magnitude compared to the non-magnetic case. It also greatly changes the dynamics 
and morphology of the protostellar accretion flow near the binary. For initially unequal mass 
binary seeds, the issue of mass ratio evolution becomes important. In §|U we follow the the 
evolution of unequal-mass protobinaries accreting from envelopes magnetized to different 
levels. We find that the well-known tendency towards equal mass in the non-magnetic case 
due to preferential mass accretion onto the secondary is weakened or even suppressed by 
magnetic braking. We summarize the main results and put them in context in § [51 



2. Problem Setup 



As mentioned in the introduction, the focus of our investigation is on the protostellar 
accretion phase of binary formation, where the magnetic effects are least explored. We 
will therefore skip the pre-stellar phase of core evolution leading up to the production of a 
pair of binary seeds, and assume that, at the beginning of our calculations, the seeds are 
already formed and are ready to accrete fro m a rotating, magnetiz ed en velope . The setup 
is similar in spirit to the influential work of iBate fc Bonnelll ( 119971 ) and iBatd (120001 ). who 
studied the effects on bina ry properties of mass accretion from an unmagnetized envelope 
(see also I Artymowica 1 1 9 831 ) . We postpone a treatment of both the pre-stellar and protostellar 
phases of binary formation to a future investigation. In what follows, we describe the initial 
conditions for the envelope and binary seeds as well as the numerical code used for following 
the envelope collapse and protobinary accretion. 

Although pro tostellar envelopes are often observed to be irregular and filamentary (e.g., 



Tobin et al. 



20101 ). it is instructive to model them simply, so that the basic effects can be 



illustrated as cleanly as possible. Since the formation and evolution of binaries involve 
complex dynamics that are challenging to simulate and interpret, it is useful to set up 
the calculations in such a way that the numerical results can be checked against analytic 
expectations. One way to achieve this is to start with a self-simila r initial configuration for 
the envelope, with an r~ 2 density distribution given by ( lShulll977l ) 

p[r) 



47rGr 2 ' v ' 

where c s is the isothermal sound speed, and A is an over-density parameter . The collapse 
of such an initial configuration is expected to remain self-similar ( 1Shulll977l ). and the self- 
similarity has proven t o provide a powerful chec k on the numerically obtained solutions (see, 
e.g., Allen et al. 2003: iMellon fc Lilbooa bood ; Kratter et al. 2010; and discussion in §EJ). 



We choose an over-density parameter A = 4 (corresponding to a ratio of thermal to 
gravitational energy of a = 3/(2 A) = 0.375), so that the initial configuration is denser than 
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the famous equilibrium singular isothermal sphere by a factor of 2. The mass enclosed within 
any radius r can be integrated as, 

Ac 2 

M(r) = -^r. (2) 

For an adopted core radius R = 10 17 cm and isothermal sound speed c s = 0.2 km/s 
(corresponding to a temperature of ~ 10 K), the above equation yields a total core mass 
Mtot = 1.2 M Q and an average free-fall time tff ~ 88 kyr. 

We generalize the isothermal configuration to include both rotation and magnetic fields. 
To preserve the self-s imilarity, the rota tion speed cannot depend on radius but can have an 
angular dependence (lAUen et al.ll2003l ). A convenient choice is = Vq sin 9 (where 9 is 
the polar angle measured from the rotation axis), which ensures that the angular speed is 
finite on the axis. We pick vq = c s , in order to have as large a rotation speed as possible, 
so that the binary can be well resolved, especially for strongly magnetized cases where the 
binary separation is reduced by a large factor compared to the non-magnetic case; even 
faster, supersonic, rotation may produce undesirable shocks. The adopted rotation profile 
corresponds to a ratio of rotational to gravita tional energy = (vn/c x ) 2 /(3A) rs 0.083, 
somewhat higher than used i n other works (e .g., 



iMachida et al.ll2010l ). Nonetheless, it is still 
within the range inferred by Goodman et al. ( 19931 ) from NH 3 observations of dense cores. 



We choose an initially uni-directional magnetic field along the rotation axis (or the z- 
axis in the simulations), with the field strength decreasing away from the axis as 1/w (where 
w is the cylindrical radius), such that the mass-to-flux ratio is constant spatially. To avoid 
singularity at the origin, we soften the profile to l/(zu + r h ), so that 



D, 



Ad 



zu 



GX^ + r h 



(3) 



where A in the denominator is the dimensionless mass-to-flux ratio of the envelope in units 
of the critical value (27rG 1 / 2 ) _1 , and the parameter is defined below. Even though dense 
cores typically have A of a few (see discussion in § [1]) , we will consider a much wider range 
of A =2, 4, 8, 16 and 32, as well as the non-magnetic case (A = oo), so as to capture any 
trend that may exist in the protobinary properties as the field strength increases gradually. 



To study the protobinary evolution, we follow iBate fc Bonnelll ( 119971 ) and iBatd (120001 ) 
and insert two equal-mass "seeds" near the center of the protostellar envelope at the be- 
ginning of the calculation; non-equal mass binaries will be treated separately in § HI We 
assume that each of the binary seeds has a small initial mass of 0.05 M Q . To determine the 
initial binary separation, we assume that the binary seeds get both their masses and orbital 
angular momentum from a sphere of radius = 8.30 x 10 15 cm in the initial envelope that 
contains 0.1 M . In other words, we assign both the mass and angular momentum inside rh 
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to the binary seeds. The initial orbital angular momentum is thus 



Th 7T 27T 



L h = L(r h ) = J J J P v o sin 9 vor 2 sin 9 drd9d(p 
ooo 

which, for circular orbits, yields an initial binary separation of 



Ac 2 s v 
3G 1 



(4) 



3.687 x 10 15 cm « 246 AU, 



9(GM b 



(5) 



where Mj, = 0.1 M is the total mass of the binary. 



To follow the e nvelope accretion and protobinary evolution, we use an MHD version 
of th e ENZO adaptive mesh refinement hydro code (IBrvan fc Norman 



(IWang fc Abel 



1997 



2009 



O'Shea et al. 



20041 ). It incorporates a sink particle treatment ( IWang et al.l l2010l ). 



The magnetic field is evolved with a conse rvative MHD solver that includes the hyperbolic 
divergence cleaning of iDedner et al.l (120021 ) . The MHD version of the code is publicly avail- 
able from the ENZO website at http://code.g oogle.com /p/enzo/\ I t has been used to f ollow 
successfully t he formation and evolution of magnetized galaxies (IWang fc Abel! l2009i ) and 
star clusters (IWang et al.l 120101 ). as well as single star formation in magnetized dense cores 
dzhao et alJboilh . 



We carry out simulations with periodic boundary conditions inside a box of length 
5 x 10 17 cm on each side, which is significantly larger than the protobinary envelope (of 
radius 10 17 cm). The region outside the envelope is filled with a uniform isothermal medium 
of the same density (and temperature) as that at the outer boundary of the envelope (p a ~ 
1.91 x 10~ 19 g cm" 3 ), so that the pressures are initially balanced across the boundary. To 
speed up computation, we adopt a relatively coarse base grid of 64 3 , but allow for 8 levels 
of refinement, with a smallest cell size of about 2.0 AU. 

As usual, we adopt a barotropic equation of state (EOS) that mimics the isothermal 
EOS at low densities and adiabatic EOS at high densities: 



P(p) = Pel 



1 



P 



Pcrit 



2/3' 



(6) 



with a critical density p 



crit 



10 



-13 g cm 3 £ Qr ^Yi e transition between the two regimes. 



We tre at the binary sta rs as sink particles. The sink particle treatment is described 

Briefly, each particle accretes according to a modified Bondi- 



in detail in IWang et al 



2010 



Hoyle formula (see iRuffert Tl994h . New sink particles are created at the center of those 
over-dense cells that violate the Jeans criterion at the highest refinement level, i.e. when 
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PrrJJ > PJ 



7.90 x 10 lA g-cm 3 , where we have used a Truelove's (ITruelove et al. 



G V NAx I 

19971 ) safety number N = 16 and a finest cell size Ax = 3.05 x 10 13 cm. Newly created sink 
particles are subject to merging, which is controlled by two parameters: the merging mass 



M, 



tii erg 



0.01 M Q and merging distance I 



merg 



5 x 10 14 cm (for details, see I Wang et al. 



20101 ). These values are chosen to eliminate artificial particles and to maximize computation 
efficiency Our main results are insensitive to these parameters as long as they are reasonably 
small. 

We have tested the sink particle accretion algorithm against known se mi- analytic solu- 
tions of the collapse of (non-magnetic) singular isothermal sphere (jShulll977l ) and found good 
agreement. For example, the collapse of our chosen initial configuration (Eq. [T| with A = 4) 
yields, after some initial adjustment, a constant mass accretion rate of 1.06 x 10~ 5 M Q yr _1 , 
which matches Shu's result almost exactly. 

In a magnetized medium, we let the sink particle accrete only the mass but not the 
magnetic field; the field is left behind. This treatment is a crude represe ntation of the 
decou pling of the mag netic field from t he matter expected at high densities (INakano et al. 
20021 ; see discussion in lZhao et al.ll201ll ). 



3. Magnetic Braking and Protobinary Orbital Evolution 



Before discussing detailed quantitative results, we first describe the qualitative effects 
of the magnetic field on the (initially equal-mass) protobinary orbit by contrasting the two 
extreme cases: the case without a magnetic field (the HD or hydro case hereafter), and 
the strongest field case (A = 2, which is consistent with available Zeeman measurements, 



see iTroland fc Crutcherl 120081 ). Both simulations reached a rather late time (t > 80 kyr, 
comparable to the free-fall time at the average density), when the binary seeds have finished 
many orbits around each other. In Fig. QQ we show the snapshots of the two cases near 
the middle point of the simulation, at t ~ 39 kyr. They look strikingly different. In the 
HD case, there are two well separated stars, each surrounded by a circumstellar disk. The 
disks are surrounded by a well-defined circumbinary structure, with two prominent spiral 
arms embedded in it. The circumstellar disks and circum binar y stru cture are similar to 



those found in previous non-magnetic calculations, such as iBatd (120001 ) and lHanawa et al. 



(120101 ). They are the result of the conservation of angular momentum in the rapid collapsing 
region of the envelope, and redistribution of angular momentum close to the binary through 
gravitational torque. 



These well-known features are completely absent in the strongly magnetized A = 2 case. 
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They are replaced by two irregular lobes, which were studied in detail in lZhao et al.l ( 1201 ll ); 
they are the so-called "magnetic decoupling-enabled structures" (or DEMS for short) pro- 
duced by the magnetic flux decoupled from the matter that enters the stars (sink particles). 
The DEMS are magnetically dominated, low-density, expanding regions. They surround 
the protobinary, whose separation is much smaller than that of the HD case (see Fig. [T|). 
Clearly, the magnetic field has greatly modified not only the circumstellar and circumbinary 
structures, but also the binary orbit. In the following subsections, we will discuss these 
modifications more quantitatively. 



3.1. Binary Separation and Angular Momentum- Mass Relation 

3.1.1. Hydro Case: Checking against Expectations 

We begin our quantitative discussion with binary separation. It is shown as a function 
of time in Fig. [2] for all initially equal- mass simulations. The evolution of the binary separa- 
tion is particularly interesting in the HD case. It decreases initially for about 20 kyr, before 
increasing almost linearly with timd^. The linear increase is a natural consequence of the 
initial protostellar envelope configuration chosen, which is self-similar over a range of radii 
(excluding the regions close to the center and the outer edge). The configuration is expected 
to collapse self-similarly, with the binary separation increasing linearly with time (for the 
same reason that the size of the expansion wave in th e well-kno wn inside out collapse of a 



singular isothermal sphere increases linear with time; IShul 119771 ). after some initial adjust- 
ment. The decrease in binary separation during the initial adjustment comes about because 
we have put all of the angular momentum of the material inside a small sphere into the 
orbit of the initial binary seeds. This leads to an overestimate of the initial binary orbital 
angular momentum (and thus the separation) because a fraction of the angular momentum 
should be left behind in the circumstellar disks and circumbinary structure. We have exper- 
imented with binary seeds formed out of a smaller sphere that have smaller initial masses 
and separation. They exhibits a similar linear growth of separation with time after a shorter 
adjustment period. In other words, the protobinary reaches the expected self-similar state 
more quickly. The agreement of the numerical result with expectations gives us confidence 
that the ENZO-based AMR code can treat the protobinary evolution problem properly. 

Another check on the numerical solution comes from the expected scaling between the 



1 The wiggles on the curve is due to small orbital eccentricity (of order 10% or less) excited during the 
protobinary evolution. We will postpone a detailed study of the magnetic effect on eccentricity to a future 
investigation. 
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Logarithm of Column Density (g/cm 2 ) Logarithm of Column Density (g/cm 2 ) 




-0.721 -0.307 0.107 0.522 0.936 1.35 1.76 -0.766 -0.525 -0.285 -0.0443 0.196 0.437 0.677 




Time= 38.882 kyr, Lengths 3.00e+16 cm 




Time= 39.457 kyr, Lengths 3.00e+16 cm 



Fig. 1. — Distribution of the logarithm of the column density E (in g-cm~ 2 ) along z-direction 
and velocity field on the equatorial plane. Both the HD (left panel) and A = 2 (right panel) 
case are at t ps 39 kyr. The well-defined circumstellar and circumbinary disks in the former 
are replaced by two magnetically dominated lobes in the latter. The binary stars are marked 
by crosses. 
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Time (kyr) 

Fig. 2. — Evolution of binary separation with time for HD (black solid), A = 32 (purple 
dotted), A = 16 (blue long-dashed), A = 8 (green dash-dot-dot-dotted), A = 4 (yellow dash- 
dotted), and A = 2 (red thick solid) cases. Note the large difference between the HD and 
the realistically magnetized (A = 2 and perhaps 4) cases. 
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orbital angular momentum and total mass of the protobinary. The scaling can be obtained 
as follows. We assume that the gas in our singular isothermal envelope accretes shell by 
shell onto the binary. The total mass M(r) and total angular momentum L(r) for a sphere 
of radius r are given, respectively, by Eq. (j2J) and 



If all of the mass and angular momentum of the material within the sphere are accreted 
onto the binary, then the orbital angular momentum and binary mass M& must be related 
through 



This relation is shown as the dashed line (scaled down by a factor of ~ 1.5) in Fig. [31 together 
with the angular momentum-mass (or Lb-Mb) relations obtained in all of our simulations. 
Again, the HD case is particularly noteworthy. Its Lb — Mb relation closely matches the 
analytical prediction from eq. ©, except for a correction factor of ~ 1.5. The correction is 
to be expected because the mass and orbital angular momentum of the protobinary at any 
given time does not come from a region of perfectly spherical shape (assumed in deriving 
the above equation). More slowly rotating material near the rotation axis can collapse more 
quickly onto the binary than that near the equator, lowering the actual angular momentum 
of the binary relative to its mass. In any case, the broad agreement between the numerical 
result on the Lb-Mb relation and the analytical expectation lends further credence to the 
correctness of the hydro simulation. 



Fig. [2] shows a general trend that the binary separation decreases with increasing mag- 
netic field strength. The difference is especially striking at late times, when the separation 
increases with time for the HD case but stays roughly constant or decreases for the mag- 
netized cases. By the time t ~ 80 kyr, the protobinary separation in the HD case reaches 
~ 500 AU, which is much larger than that in the A = 4 (~ 50 AU) and A = 2 (~ 10 AU) 
case. There is little doubt that a realistic magnetic field (corresponding to a dimensionless 
mass-to-flux ratio of a few, see Troland & Crutcher 2008 and discussion in § [T|) can shrink 
the protobinary orbit by a large factor. 

The reduction in binary separation is related to a decrease in the orbital angular mo- 
mentum of the system. This is shown explicitly in Fig. [3j There is a general trend for the 
orbital angular momentum to decrease with increasing magnetic field strength. When the 
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3.1.2. Magnetic Effect on Protobinary Orbit 
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Fig. 3. — The relation between the orbital angular momentum and total mass of the proto- 
binary. The dash line is the the analytical prediction from Eq. [8] (reduced by a correction 
factor of ~ 1.5). The different curves are: HD (black solid), A = 32 (purple dotted), A = 16 
(blue long-dashed), A = 8 (green dash-dot-dot-dotted), A = 4 (yellow dash-dotted), and 
A = 2 (red thick solid). 
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mass of the stellar seeds quadruples from 0.1 to ~0.4 M , the orbit angular momentum in 
the HD case increases by more than an order of magnitude, whereas that in the strongly 
magnetized A = 2 case increases by merely a factor of 2. In the latter case, there is a large 
amount of mass added to the protobinary but relatively little angular momentum. Since the 
binary separation a is sensitive to the orbital angular momentum Lb (i.e., a oc L\ for a fixed 
binary mass M&), even a relatively modest change in the angular momentum would lead to 
a significant change in the separation. 

We conclude from Figs. [3] and |5] that the main magnetic effect on the protobinary 
evolution is to reduce its orbital angular momentum compared to the non-magnetic case, 
which in turn leads to a tighter orbit. This result may appear puzzling at the first sight, 
because the magnetic forces do not act directly on the binary seeds. However, they do act 
on the material to be accreted by the seeds. By changing the angular momentum of such 
material, the magnetic field can greatly affect, perhaps even control, the orbital evolution of 
the protobinary. In the next subsection, we explore in some detail the mechanism through 
which the magnetic field shrinks the protobinary orbit. 

3.2. Magnetic Braking and Angular Momentum Removal from Protobinary 



It is well-known that magnetic fields interact strongly with fluid rotation, through mag- 
netic braking. If the infalling material has a large fraction of its angular momentum removed 
prior to its arrival at the binary seeds, it would add mass but relatively little angular mo- 
mentum to the protobinary system. As a result, the binary separation would increase less 
rapidly with time compared to the hydro case; it may even decrease with time if the magnetic 
braking is strong enough. 

The presence of magnetic braking can be seen directly from Fig. HI where we plot the 
field lines at the representative time t ~ 39 kyr for the A = 2 case. The initially straight field 
lines are twisted by the fluid rotation into helical shape, especially in the polar regions. The 
magnetic tension force associated with the twist acts back on the fluid, braking its rotation. 

To quantify the strength of the magnetic braking, let us consider a finite volume V with 
surface S. Inside this volume, the total magnetic torque relative to the origin (from which a 
radius vector r is defined) is 



Accretion Flow 




(9) 



Fig. 4. — 3D view of the magnetic field lines and an iso-density surface in the inner part 
of the protobinary accretion flow at the representative time t ~ 39 kyr, showing the field 
twisting that is the smoking gun of magnetic braking. The plotted region has a dimension 
of 5 x 10 15 cm. 
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where the integration is over the volume V. Typically, the magnetic torque comes mainly 
from the magnetic tension rather than pressure force. The dominant magnetic tension term 



can be simplified to a surface integral (IMatsumoto fc Tomisakal 120041 ) 



N * = 47 /( rxB X B - rfS )> ( 10 ) 

over the surface S of the volume. This volume-integrated magnetic torque is to be compared 
with the rate of angular momentum advected into the volume through fluid motion, 

N - = ~J P(rx v)(vrfS), (11) 

which will be referred to as the advective torque below. 

Since the initial angular momentum of the protobinary envelope is along the z-axis, we 
will be mainly concerned with the z— component of the magnetic and advective torque, 



N t , z 

and 



^ J(xB y -yB x )(B-dS) (12) 



N a , z = - J p{xv y - yv x )(v ■ dS), (13) 
which can change the z— component of the angular momentum L z within the volume V. 

As an example, we show in Fig. |5] the distributions of the magnetic and advective 
torques N t<z and N a>z for cubic boxes of different sizes that are centered at the origin, at 
the representative time t ~ 39 kyr for the A = 2 case. As expected, the volume-integrated 
magnetic torque is negative for boxes of most sizes; it removes angular momentum from the 
material inside the boxes through magnetic braking. At the time shown, the braking torque 
is particularly large outside ~ 10 16 cm from the origin, reaching absolute values of order 
10 41 g cm 2 s -2 or larger. To appreciate how large this torque is, we note that the orbital 
angular momentum of the protobinary in the hydro case is ~ 4.4 x 10 52 g cm 2 s _1 around the 
same time. This angular momentum would be removed by the above magnetic torque on a 
time scale of 13 kyr, much shorter than the envelope collapse time. Given the magnitude of 
the magnetic torque, it is not surprising that the orbital angular momentum (and thus the 
separation) of the protobinary can be reduced significantly. 

The magnetic torque changes with time, however. At the time shown in Fig. the 
magnetic torque overwhelms the advective torque for boxes of most sizes, leading to a net 
decrease of the angular momentum with time for the material in these boxes. This may not 
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Time= 39.457 kyr 




1Q 15 1Q 16 1Q 

Box Half — width (cm) 



Fig. 5. — The magnetic (black '+') and advective (blue 'x') torque and the sum of the two 
(red) for cubic boxes of different half-width b for the A = 2 case, at a representative time 
t pa 39 kyr. A positive torque increases the angular momentum within a volume whereas a 
negative one decreases it. 
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happen at other times. To evaluate the accumulative effects of the magnetic and advective 
torque over time, we define for each box of half width b 

L tjZ (b,t) = [ N tjZ (b,t')dt' (14) 
Jo 

and t 

L a , z (b,t) = [ N ajZ (b,t')dt' (15) 
J o 

which are the amount of the z— component of the angular momentum inside the box changed 
by the magnetic and convective torque, respectively, up to time t. These two quantities are 
to be compared with the actual change of the angular momentum inside each box between 
time t and t = 0, 

AL z {b, t) = L z {b, t) - L z {b, t = 0). (16) 

Barring a significant gravitational torque, one expects AL z (b,t) to be close to L ajZ {b,i) 
in the hydro case, because fluid advection should be the main channel for angular momentum 
change. This is indeed the illustrated in panel (a) of Fig. [61 where both AL z (b, t) and 

L a ,z{b, t) are plotted as a function of the box size b at the representative time t « 39 kyr. The 
change in angular momentum AL z (b, t) does follow closely the advected angular momentum 
for boxes of sizes larger than ~ 10 16 cm. For boxes of smaller sizes, there is significantly 
more angular momentum advected into a box than the actual angular momentum change 
in it, indicating that a good fraction of the angular momentum advected into the box is 
transported back out, presumably by the gravitational torques associated with the spiral 
arms that are visible in the left panel of Fig. [U The gravitational torques may also be 
responsible for the small excess of L a ^ z (b,t) over AL z (b,t) for the larger boxes. 

The situation is very different in the presence of a relatively strong magnetic field, as 
shown in panel (b) of Fig. El Plotted are L t , z (b, t), L a>z (b, t), and L tiZ (b, t)+L ajZ (b, t) together 
with AL z (b, t) as a function of box size at t ~ 39 kyr for the A = 2 case. The total change of 
the angular momentum due to the magnetic and advective torque acting on the boundary 
of a box, L tz (b, t) + L a>z (b, t), is very close to the actual change in angular momentum inside 
the box, AL z (b,t), indicating that any additional torques (such as gravitational torques), if 
present, play a relatively minor role in angular momentum transport. This is not surprising, 
because the prominent spiral arms of the hydro case are disrupted by the magnetic field 
completely. More importantly, the sum L t<z (b,t) + L a . z (b, t) has a magnitude much smaller 
than \L ttZ (b, t)\ and \L ajZ (b, t)\ individually, which means that most of the angular momentum 
advected into a box is removed by the magnetic braking, leaving little net angular momentum 
change in it. This is unequivocal evidence that the magnetic braking in the A = 2 case 
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Fig. 6. — The integrated magnetic and advective torques L tjZ (black solid), L a)Z (blue dotted), 
and their sum L t>z + L a>z (red dashed) together with the actual angular momentum change 
AL (green thick solid) for cubic boxes of different half-width b. The HD case is shown in 
panel (a) with zero L t , z and the A = 2 case is shown in panel (b), both at a similar time 
t ~ 39 kyr. The vertical dashed line in panel (a) indicates the approximate position of the 
sink particles, whereas in panel (b) such particle position lies below 10 15 cm. 
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controls the angular momentum evolution of the protostellar accretion flow, which in turn 
shapes the orbit of the protobinary. 

The effect of magnetic braking is even stronger than indicated by panel (b) of Fig. O 
This is because L a)Z (b,t) is the net angular momentum advected into a box, i.e., the sum 
of the positive angular momentum advected into the box L+ Z (b,t) and the negative angular 
momentum advected out of the box L~ z (b, t). In the A = 2 case, the magnetic braking drives 
an outflow, which produces a negative angular momentum L~ z (b, t) that is not much smaller 
in magnitude than L+ Z (b,t). In other words, the angular momentum advected by infall into 
a box is larger than that shown in Fig. [61 and most of this larger angular momentum is 
removed by both the magnetic braking itself and the braking-induced outflow. 

To illustrate the effect of magnetic braking more visually, we plot in Fig. [7] the distribu- 
tion of specific angular momentum on the equatorial plane for the hydro and A = 2 case at 
the representative time t ~ 39 kyr. It is clear that, for the hydro case, the specific angular 
momentum is roughly constant within a dimensionless radius of ~ 0.05 (or ~ 2.5 x 10 16 cm), 
indicating that the collapsing material has a more or less conserved angular momentum be- 
fore it is accreted by the protobinary. The relatively large specific angular momentum of 
the accreted material is what drives the binary separation to increase. In the A = 2 case, 
the specific angular momentum of the material to be accreted is much smaller; it is reduced 
by twisted field lines (see Fig. H|) as the material falls toward the binary. It is the accretion 
of the severely braked, low angular momentum material that drives the protobinary closer 
with time. 



3.3. Circumstellar and Circumbinary Structures 

We have already seen from Fig. [1] that a realistic magnetic field of A = 2 changes the 
circumstellar and circumbinary structures of the HD case greatly: the prominent features of 
the HD case, two well-defined circumbinary disks and a circumbinary disk with two promi- 
nent spiral arms, are replaced by two low density lobes, which are filled with the magnetic 
flux decoupled fr om the matter that has been accreted onto t he binary, i.e., the DEM S 



(IZhao et al.ll201ll ; see also Seifried et al. 2011; Joos et al. 2012; iKrasnopolsky et al.ll2012l ). 
These strongly magnetized structures present an obstacle to mass accretion onto the pro- 
tobinary. Although the DEMS tend to be less prominent for weaker magnetic fields, they 
still dominate the face-on view of the A = 4 and 8 cases, especially at late times. An im- 
plication is that, for protostellar envelopes magnetized to a realistic level (with A ~ a few), 
the protobinary mass accretion does not follow the traditional path: from the envelope to a 
circumbinary disk to circumstellar disks to individual stars. Rather, the envelope material 
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Fig. 7. — Snapshot of the gas specific angular momentum (in g ■ cm 2 ■ s _1 with logarithm 
scale) on the equatorial plane of the inner accretion flow for the HD case (left panel) and 
the A = 2 case (right panel) at t ~ 39 kyr, showing the strong braking of the material to be 
accreted by the protobinary in the magnetized case compared to the HD case. The arrows 
are velocity vectors, and the crosses mark stars. Only the central region of 10 17 cm on each 
side is plotted. 
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collapses directly close to the stars typically, and be accreted along azimuthal directions not 
occupied by the DEMS. 

Another magnetic effect on the circumbinary environment is illustrated in Fig. [HI where 
we display both the edge-on and face-on view of the A = 2 case at a relatively early time 
t = 18 kyr. The edge-on view shows that the protobinary is surrounded by two large, mostly 
expanding, regions (one each above and below the equator) that are absent in the HD case. 
Their dynamics are dominated by the rotationally twisted magnetic fields (see Fig. H]). It 
is in these regions that most of the angular momentum extracted magnetically from the 
material that falls into the stars is stored. Such bipolar expan ding; regions are seen in many 



magnetized core collapse simulations (for an early example, see lTomisakalll998l ) . They block 



the protobinary accretion over most of the solid angle, and force the accret ion to occur 



main ly through a flattened, equatorial structure — a circumbinary pseudodisk (jGalli fc Shu 



1993|). 



The pseudodisk is also clearly visible in the face-on view of the system in the right panel 
of Fig. [HI as a nearly circular region of enhanced column density. Unlike the circumbinary disk 
in the HD case (see the left panel of Fig. [T]), the pseudodisk is not rotationally supported. 
It falls supersonically inward, and the infall can be seen in the velocity fields shown in 
Fig. [HI Another difference is that the circumbinary pseudodisk does not have prominent 
spiral arms. This is not surprising, because the pseudodisk is strongly magnetized, which 
makes gas compression more difficult. The rapid infall also leaves little time for the spirals 
to develop. The lack of spiral arms is consistent with our earlier result that the angular 
momentum transport in the A = 2 case is dominated by magnetic braking rather than 
gravitational torque. 



3.4. Protobinary Mass Growth 

Despite the differences in the morphology and dynamics of the circumbinary environ- 
ment with and without a magnetic field, the stars grow at remarkably similar rates in all 
cases (of order ~ 10~ 6 M Q /yr). This is shown in Fig. [9], where the mass of each component 
of the protobinary is plotted as a function of time, for different degrees of magnetization. 
The reason for the similarity is that the mass accretion rate is determined mainly by the 
dynamics on the envelope scale: the material collapsing down from the envelope will find a 
way to the stars sooner or later, irrespective of the details of the circumbinary environment. 
In the HD or weakly magnetized cases, there is more material parked in the circumstellar 
and circumbinary disks, making the mass accreting rate onto the stars lower. However, this 
effect is smaller than that from the retardation of envelope collapse by magnetic tension in 
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Fig. 8. — Left panel: edge-on view of the column density in the A = 2 case at t — 18 kyr, 
showing two slowly expanding polar regions sandwiching an infalling, circumbinary pseu- 
dodisk in the equatorial region. Right panel: face-on view of the same structure, showing 
that the circumbinary pseudodisk is not rotationally supported; it collapses rapidly. 
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the stronger field cases. As a result, the protostellar mass accretion rate is somewhat lower 
in the stronger field cases, although not by a large factor. 

It is reassuring that the masses of the two initially equal-mass components stay nearly 
the same for all cases, despite the fact that the inner protobinary accretion flow can become 
rather disordered at times, especially for strongly magnetized cases. The magnetic field does 
not appear to change much the rate of mass accretion by one component relative to the 
other, although the situation is drastically different for initially unequal mass protobinaries, 
as we demonstrate next. 



4. Magnetic Braking and Mass Ratio of Unequal-Mass Protobinaries 



The vast majority of low (Sun-like) mass binaries are unequal- mass systems (jDuquennoy fc Mayor 



199ll ; iRaghavan et al.ll2010l ). The mass ratio, defined as q = M 2 /Mi where Mi and M 2 are 
the mass of the primary and secondary respectively, is one of the fundamental parameters 
that characterize the binaries. Its observed distribution provides important constraints on 
binary formation and evolution (see discussionn in §EJ)- 

The mass ratio q of protobinaries is expected to be affected by magnetic braking. This 
is because the change of q depends on the specific angular momentum of the circumbinary 
material to be accreted relative to that of the binary. It is well known that high angular 
momentum material tends to accrete preferentially onto the less-massive secondary, which 
has a higher specific angula r momentum than the primary, driving the mass ratio toward 



unity (IBate &: Bonnelll 119971 ; iBate et al.ll2002l ). As we have seen above, magnetic braking 



can remove the angular momentum of the circumbinary material efficiently. It is expected 
to weaken the tendency for preferential accretion onto the secondary. We show that it is 
indeed the case in Fig. [TOj 

Shown in Fig. [TU]are the mass ratios as a function of time for three cases A = oo (HD), 
16 and 2. The initial conditions are the same as those for the equal-mass cases discussed 
in § El except that the mass and angular momentum inside the radius of the power-law 
envelope are given to a binary system of mass ratio q = 0.25 rather than q = 1. The mass 
ratio increases with time in the HD case, in agreement with previous results. As expected, 
the magnetic field slows down the increase in q compared to the HD case. Even a rather 
weak magnetic field of A = 16 has an appreciable effect on the mass ratio evolution. For 
the more realistic value of A = 2, the primary accretes much faster than the secondary, 
by a factor of 3-4, so that the mass ratio remains roughly constant. The large contrast 
between the HD and A = 2 case leaves little doubt that magnetic braking is an important 
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Fig. 9. — Stellar mass (in solar units) growth of the initially equal-mass binary system for 
HD (black solid), A = 32 (purple dotted), A = 16 (blue long-dashed), A = 8 (green dash-dot- 
dot-dotted), A = 4 (yellow dash-dotted), and A = 2 (red thick solid) cases. Note that, in all 
cases, the two stars stay roughly equal-mass at all times. 
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Fig. 10. — Evolution of the protobinary mass ratio q with time for envelopes that have 
different levels of magnetization, with A = oo (HD, black solid), 16 (blue dashed), and 2 (red 
thick solid). 
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factor to consider in understanding the mass ratio distribution of binaries. We will leave a 
comprehensive exploration of this issue to a future investigation. 



5. Summary and Discussion 

We have carried out a set of idealized numerical experiments to demonstrate the effects 
of the magnetic field on protobinary evolution during the mass accretion phase. The pro- 
tostellar envelope was idealized as a rotating, magnetized singular isothermal sphere. Its 
collapse onto a pre-existing pair of binary seeds was followed using an MHD version of the 
ENZO AMR hydro code that includes a sink particle treatment. We found that a magnetic 
field of the observed strength (corresponding to a dimensionless mass-to-flux ratio of a few) 
can remove, through magnetic braking, most of the angular momentum of the material that 
reaches the protobinary. The reduction in the angular momentum of the protobinary ac- 
cretion flow has two important consequences: compared to the non-magnetic case, (1) the 
protobinary orbit becomes much tighter, and (2) the mass-ratio does not increase as fast with 
time for initially unequal mass systems. In addition, the magnetic field drastically changes 
the morphology and dynamics of the structures that surround the protobinary. It suppresses 
the formation of the familiar circumstellar and circumbinary disks in the non-magnetic case, 
as well as the spiral arms embedded in them. These structures are replaced by a bipolar 
magnetic-braking driven expanding regions, a dense, infalling, circumbinary pseudodisk in 
the equatorial region, and low-density, highly magnetized structures close to the protobinary 
that expand against the pseudodisk (the DEMS). We conclude that both the basic charac- 
teristics (such as mass ratio and separation) of the protobinaries and the environment in 
which binaries grow are strongly modified by a realistic magnetic field. 

The magnetic braking-driven inward migration of protobinaries may have observable 
consequences. As mentioned in the introduction, the distribution of binary orbital separation 
during the earliest, Class phase of star formation may be different from those at later times. 
In particular, there is tent ative evidence for a "desert" free of Cla ss binaries with separation 



between ~ 150 — 550 AU (IMaury et al.ll2010l ; Enoch et al.ll201ll ). which is not present in the 



Class I or later phases. If confirmed by future observations, this gap must be filled in by 
binaries from either outside the gap or interior to it. The magnetic braking is an efficient 
way of shrinking the protobinary orbit. It can in principle move some binaries born on wide 
orbits outside the gap into the gap. Indeed, Offner et al. (2009, 2010) found that most of the 
binaries in their cluster formation simulations were born with relatively wide separations, as 
a result of radiative feedback, which tends to suppress close binary formation through disk 
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fragmentation; the orbits of such wide binaries could be tigthened by magnetic brakingp. 
A potential problem is that the braking may be so efficient during the protobinary mass 
accretion phase (when a massive, slowly rotating envelope is still present) that the binary 
separation would move quickly through the gap and pile up below the resolution of the 
current generation of millimeter/submillimeter interferometers (~ 50 — 100 AU). If this is 
the case, one would expect to find an over- abundance of relatively close protobinaries with 
separation < 100 AU that may be uncovered by ALMA and JVLA. The same pileup at small 
separations would also occur if the binaries born on tight orbits interior to the gap are kept 
from expanding into the gap by magnetic braking. 

The problem of potentially turning most wide binaries into close binaries is related to 
the magnetic braking catastrophe in disk formation (Galli et al. 2006). Both semi-analytic 
arguments and numerical simulations have shown that, in magnetized laminar dense cores 



braking in the ideal MHD limit ( 


Allen et al. 


2003 




Galli et a 


2006 




Price & Bate 


2007 


Mellon & Li 


2008 




Hennebelle & Fromangl 


200 


8; 


Seifried et al. 


2011 




Dapp & Basu 


2011) 



The magnetic braking must be weakened somehow in order to form both large-scale disks 
and wide binaries. In the context of disk formation, there are several proposed mechanisms 
for weakening the magnetic bra king, includin g the misalignment between the magnet i c and 



rotation axes (jJqos et al.l |2012| ) , turbulence fjSantos-Lima et al.l 120121 ; ISeifried et al.l 12012 



Myers et al.ll2012l ). and the depletion of the slowly rotating pr otostellar env elope that acts 



to br ake the disk, either by outflow stripping ( jMellon fc Lill2008l ) or accretion ( iMachida et al. 



2010|). The last possibility is particularly intriguing, because it implies a rapid growth of the 
rotationally supported disk during the transition from the deeply embedded Class phase to 
the more revealed Class I phase that can be observationally tested. Similarly, the depletion 
of the protobinary envelope and the associated weakening of magnetic braking may enable 
the separation of the protobinaries to grow quickly during the Class 0-Class I transition. 
If this is the case, the orbits of wide binaries may first shrink during the Class phase 
due to efficient magnetic braking and then re-expand as the protobinary envelope depletes. 
This and other possibilities for protobinary migration should be testable with high-resolution 
millimeter/submillimeter interferometric observations, especially using ALMA and JVLA. 

Our calculations demonstrated that magnetic braking i s important for t he evo lution of 
not only the binary separation, but also the mass ratio q. iRaghavan et al.l (120101 ) found a 



2 Batc (2012) found that the binaries in his radiation hydro simulations of cluster formation have properties 
consistent with observations, indicating that magnetic effects are not needed. It would be interesting to 
quantify how a magnetic field of the observed strength modifies the properties of the binaries formed in a 
cluster environment. 
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roughly flat distribution for q between 0.2-0.95, with a valley below 0.2, and three spikes 
around 0.25, 0.35 and 1.0 (see the left panel of their Fig. 16). It is tempting to attribute the 
last spike at q ~ 1 to the preferential accretion of high specific angular m omentum mate- 



rial onto the lower-mas s component found in hydrodynamical simulations ( iBate fc Bonnell 



1997 



m 



Bate et all 120021) . However, the spike around q ~ 1 is just o ne of the three spikes 



Raghavan et al.f s data, and it does not show up in the sample of iDuquennoy &: Mayor 
( 1199ll ) at all (see their Fig. 10). In any case, there are many more systems with mass ratio q 
well below unity than close to unity. A mechanism must be found to prevent the majority of 
the low-mass ratio protobinaries from becoming equal-mass systems. The magnetic braking 
highlighted in this paper is one such mechanism. It is a key factor to consider in under- 
standing binary formation and evolution in dense cores that are observed to be significantly 
magnetized. 

We thank P. Arras, S. Offner, K. Kratter, C. Matzner, and A. Maury for useful discussion 
and P. Wang for advice on the ENZO code. This work is supported in part by NASA 
NNX10AH30G. 
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